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O ■ Abstract: The definition of the pairing interaction for rotating systems is 

discussed, and a simple, mainly analytic model for its treatment is formu- 
lated. The results are used to describe microscopically the giant angle dipole 
resonance of the superfluid deformed nuclei in the semiclassical approxima- 
tion of the cranking model. 
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I. Introduction 



CN 
> 
CN 

in 
o 

^ . The microscopical description of the superfluid rotating nucleus has been 

CO ■ the object of a constant interest during the time, started after the early at- 

tempts, made in solid state physics, to treat the similar situation of supercon- 
ductors placed in magnetic fields [1] [2]. This similarity is better shown in the 
semiclassical approximation given by the cranking model, where the angular 
velocity from the Coriolis term is the correspondent of the external magnetic 
field. For the deformed nuclei, initially the main problem was to obtain the 
quasiparticles of a hamiltonian consisting of a cranked single-particle Nils- 
son term, with pairing interaction [3]. The further studies were concentrated 
on a more elaborated treatment of the symmetries, using the Hartree-Fock- 
Bogolyubov (HFB) solutions of the cranking hamiltonian in order to achieve 
new trial wave functions by angular momentum [4] or particle number [5] 
[6] projection. Considering the cranked HFB theory as a large-deformation 
limit, for a more accurate treatment by variation after projection [7] [8], the 
spurious rotational energies were also studied [9]. 

Beside the problems raised by the breaking of the symmetries, there is 
also another important point, related to the choice of the single-particle ba- 
sis used to define the pairing interaction term for systems which do not have 
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a time-reversal invariant single-particle hamiltonian. As it was remarked 
before in the case of the electrons in solids with non-magnetic impurities, 
the correct choice is given by the eigenfunctions of the whole single-particle 
part of the hamiltonian. This part includes the perturbation terms, and the 
final result cannot be reached simply by re-writing the pairing term from 
the unperturbed basis to the new one [1]. For nuclei a similar treatment 
was given in [4], where the pairing interaction defined in the Nilsson basis 
is re-written numerically in the basis where the Nilsson plus Coriolis single- 
particle term is diagonal. Later on, a different choice was presented in ref. 
[5], and used to estimate the increase in the moment of inertia due to the 
Coriolis anti-pairing and rotational alignment phenomena. At backbending, 
the quasiparticle vacuum was proved to become unstable, with effects simi- 
lar to the gapless superconductivity already known from the reference [10], 
where the Nilsson basis was used. 

In this paper the basis problem is treated in a way similar to the one 
presented in [4], but with two exceptions. First, the pairing interaction is 
defined in the new basis by using the time-reversal operator, and second 
the single particle hamiltonian is restricted to the cranked anisotropic har- 
monic oscillator (CO.) term. This term can be diagonalized analytically [11] 
[12], and moreover, the restriction is able to give a model where the relevant 
physics of the spin-independent phenomena combining the effects of defor- 
mation, pairing and rotation is maintained in the simplest way. 

The interest in such a treatment is twofold, because it gives us a better 
understanding of the pairing interaction, and also because the flow pattern, 
easy to obtain by a cranking calculation [13]- [15], is important for the mi- 
croscopical description of the isovector angle vibrations (A.V.) [16]. Even 
though the picture of two ellipsoids in opposite angular rotational oscilla- 
tion, initially proposed for these modes [17], was subject to criticism [18] 
[19], it seems to be the appropriate classical analog for the A.V. states [20] 
[21]. This picture is derived also in the classical limit of the IBA-2 model 
[22]. As it will be shown in Sect. Ill, if the contra-rotation is considered 
slow with respect to the intrinsic motion of the nucleons, and the resulting 
time-dependent wave function is properly quantized, then a semiclassical de- 
scription of the A.V. modes in the framework of the cranking model becomes 
possible. 

Along the lines sketched above, Sect. II contains the treatment of the 
pairing interaction by using the HFB trial wave functions as SO(2n) co- 
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herent states. As a counterpart, in Appendix the sp(4,R) structure of the 
single-particle hamiltonian is emphasized, and closed formulas for the matrix 
elements of the time-reversal operator are given. This algebraic formulation 
is intended to allow a direct link between the present A.V. cranking model 
and the models based on boson expansions as IBA-2. 

In this framework one can also do a further accurate treatment of the 
pairing rotation and vibration modes in the rotating nuclei by global group 
theoretical methods, giving new insight on the recently predicted pair trans- 
fer [23]. 

The details of the proposed formalism describing semiclassically the A.V. 
states in realistic nuclei will be presented in Sect. III. Finally, in Sect. IV 
the parameters are fixed, and some numerical results concerning the ground- 
state rotational band and the A.V. modes in rare-earths nuclei are discussed. 

II. The Ground State Band 

In the CO. approximation for the single-particle part of the hamiltonian, 
the writing of the pair operator pt is difficult because the basis eigenfunctions 
^nin 2 n 3 (eq. A. 8) are not time- reversal (complex conjugation) invariant due 
to the cranking term, and their complex conjugates ^* nin2n ,, are not eigen- 
functions due to the anisotropy. In this case the usual pair operator P\ 
constructed by the pairing of the time-reversed states, will be replaced by a 
more general one, given by: 

Pt = £</J|£ |a)ct cj = £ s a q ba clcl a , (1) 

a/3 a,b,a 

where a = (ni, n 2 , n 3 ), a = (a, a), \a) = \E r ni n2n 3 | cr )> a — ±1/2 is the z- 
axis spin component, s a = (— l)^~ a , K is the time- reversal operator, and 
q ba = (^ff,\K\'^f a ) = This expression originates in the commonly 

accepted form of a separable interaction term in the particle-particle (pp) 
channel, H int = ~xP'P- Here P^ is written in a basis independent form 
by using a single-particle operator Q as P^ = J2 a ,/3(P\Q\ a ) c a c 8 ( re f- [24] p. 
180). If Q = Kr l Yi m then the usual form of the multipole pairing operator is 
recovered, and in particular, when 1 = 0, this reduces up to a constant factor 
to(l). 

The ground state wave function \Z) for the proton or neutron system 
is variationally obtainable, looking for the minimum of the energy E = 
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(Z\H\Z)/(Z\Z) with the constraint 



<*> = iffi = ^ (2) 



imposed by a fixed mean-value N part of the particle number operator iV 



J2aa c l,a c aa- The many-body hamiltonian H is defined by 



H = J2e a clc aa - ^Ptp (3) 



where G is the pairing constant and e a are dependent on the fij frequencies 
(A.6,7) through 

3 1 

Criira 2 n3 = + — ) • (4) 

If uj = the matrix g = [g ab ] is diagonal, and the state |Z)| W=0 has a BCS 
form: 

|BCS> = n a ^=4=e^ p V> , ^i = E^Lc„- CT . (5) 

V 1 + Z a Z a a 

When u > the degeneracy of the energies e a changes, the matrix q is no 
more diagonal, and HFB trial functions must be used. For simplicity it is 
useful to treat these functions as elements of the Kahler manifold [25] of 
coherent states for the SO(2n) group: 

\Z) = ^^e^^ z ^\0 p ) (6) 

C = i + zz^ , z ] = (z T y 

parameterized by the antisymmetric complex matrix Z = [z a p\, and hav- 
ing the bare vacuum |0 P ) as the highest weight vector [26] [27]. The in- 
teger n specifying SO(2n) represents the number of single-particle states 
a = (ni,n2,ri3,a) under consideration. In this formulation the calculus 
of the mean values {Z\d a c a \Z), (Z\cl,Cp\Z) and {Z\(^ a CpC^cs\Z) from the 
energy function E is reduced to some derivatives of y/det( with respect 
to the coordinates z a p. Denoting by A = G{Z\P^\Z) /2, e a p = e a S a bS aa i, 
Qp a = s a q ba 5 a >- a , a = (a, a), (3 = (b,a'), Tr(A) = £ Q A* Q , the expectation 
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values (Z\Y,ap e ap^ a cp\Z) , (Z\P^\Z) and G{Z\P^P\Z) /A are, respectively 
Tr[e(l - C _1 ) T ], Tr^C^Z] and 

^ + ^[qq* - QrV - qcV^c 1 ^] • (7) 

If in the last formula only the first term | A| 2 /2 is retained, and the constraint 
(2) is accounted by introducing the Fermi level A, then the variational equa- 
tions d(E — \{N)) / dz* a p = can be explicitly written in a matrix form as 

A*ZQZ - (fZ + Zf) + AQ* = , (8) 

where / = e — Al. This equation can be further simplified because the matrix 
/ is a constant in the spin space, and Q factorizes as q <g> s, with the matrix 
Q — [Qab] real, symmetric, and s aa i = {—l)^~ a 8 a i- a . Thus, by choosing Z 
as a direct product Z = t <g> s (reduction from SO(2n) to USp(n) coherent 
states [28] ) with t real and symmetric, and denoting by trA = ^2 a A aa , the 
equations (8), (2) become 

iqi+^(ft + tf)=q (9) 

2tr[P(l + P)- 1 ] = N part . (10) 

Once we know the t matrix the ground state problem is solved, but the de- 
scription of the quasiparticle excitations, based on the calculus of the matrices 
U and V from the HFB transformation 

cL = E UpaC^ + VfaCp (11) 
3 

is not completely achieved. The correspondence between the HFB and co- 
herent states formulations becomes more evident by noticing that (8) follows 
from the defining equation of U and V at the stationary point [24] p. 613, 
[29] : 

A*U j QU* - {U ] fV* - V ] f T U*) - AV^Q*V* = (12) 

if U is non singular and Z = VU -1 . 

The system (9), (10) can be solved separately for protons and neutrons, 
obtaining anw- dependent matrix t. To construct the intrinsic states of the 
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ground state rotational band by using the functions \Z p (u p )), \Z n (u n )) deter- 
mined by the corresponding matrices t, we suppose the interaction between 
protons and neutrons being strong enough to prevent any relative motion of 
their deformed mean-fields. Consequently, the angular frequencies u p and uj n 
are equal. For a state with a given spin J, their common value ujj is obtained 
from the semiclassical condition ( [24] p. 130) : 

(Zp\L xp \Z p ) + (Z n \L xn \Z n ) = \J 1(1 + 1) ,7 = 0,2,4,... (13) 

imposed on the expectation values of the many-body angular momentum 
operators L xp , L xn . The rotational energy in the laboratory frame is the sum 
between the proton and neutron terms, 

Ei = SpM +£ n (oJi) (14) 

both defined by: £ T {u) = E T (uo) - E T (0), 

E T (u) = {Z T \H\Z T )( u )+u{Z T \L XT \Z T )( u ) ,r = p,n. (15) 

As a function of u>, E T will be interpolated in applications by using the 
following relation 

E T (u) = E 0T + A T \u\ + 1 -^ . (16) 

Here Eq t represents the constant intrinsic energy and I T is interpreted as 
the dynamical moment of inertia associated with the classical free rotation 
around the x-axis. 

III. The Angle Vibrations 

For an elastic force acting between the proton and the neutron ellipsoids 
[17], at the classical level a potential energy term 2 C$ 2 /2 depending on the 
angle $ formed by the symmetry (z) axes of the ellipsoids must be added 
to E — E p (ujp) + E n (uj n ). Considering only rotations around the x-axis, the 
Lagrange function is, up to a total time- derivative, given by 

IpUJ 2 p I n u 2 n C$ 2 . . 

£ = + — g ,u p = (f)p ,u n = (f) n . (17) 

2 this term can be also derived by expanding in $ the mean value — x{9\Qp'Qn\g) of the 
microscopic quadrupole p-n interaction [34] taken on the laboratory functions \g), (23). 
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If new coordinates {$, 0} are introduced by the transformation 



$ = 0p"0n , 0=^L±Mi , (18) 

ip ~r in 

then C can be written as a sum between a rotational and a vibrational part: 

i + 1 ■ ii <i> 2 c$ 2 

c= h±^^ + Jv±n_^L_<^L . (19) 

2 _Zp ~\~ In 2 2 

In these coordinates the classical Lagrange equations have a simple solution 
for the motion of the proton and neutron mean-fields, representing an angu- 



lar vibration with the frequency fl = yo(/ p + I n )/I p I n superposed over an 
uniform rotation with the angular velocity uo r : 

4> p (t) = u r t + a p sin fit , (fin(t) = uj r t — a n sin fit , (20) 

ujp{t) = oj r + A v cos fit , u) n (t) = uj r — A n cos fit , (21) 

where 

A p = fla p = " uj , A n = fla n = uj . (22) 

ip ' i-n ip i in 

If the angle vibration takes place adiabatically, recalling the results of the 
cranking calculations presented in Sect. II, the wave functions of the ellipsoid 
t = p,n, in the intrinsic system and in the laboratory frame, are \Z T ),^ and 
exp(— i(p T L XT )\Z T ) respectively. Considering only the case u r = 0, for 
the whole nucleus a periodic time- dependent wave function \g)(t) can now be 
obtained: 

\g)(t) = e-^-^ L -\Z p ) M \Z n ) M . (23) 

This function, subject to a method of quantization, will be used to give an 
approximation to the A.V. stationary state. The quantization problem for 
the time-dependent solutions given by the various semiclassical methods em- 
ployed in the microscopic models [30] [31] and also for some of the classical 
dynamics relevant in the description of the collective modes in nuclei [32] [33] 
was much discussed, even if an unified picture is not yet established. As it 
was remarked in [31], the wave functions generated via the time-dependent 
variational principle from trial functions with fixed norm are determined only 
up to a time-dependent phase factor. Supposing this to be true also for |^)( t ) 
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given above, the factor will be found by requiring an optimal time-evolution 
of the function \ip)(t) = a (t) ete ^\g)(t) relatively to the total microscopic hamil- 
tonian H [34]. More precisely, the integral 

X^ = J dt(^\id t - H|V> (24) 

must be stationary at infinitesimal variations of the functions a and 6, and 

5 a ,e^ = . (25) 

A first step towards finding explicitly the stationary wave function is made 
by retaining from the solution \ip) of this equation only the part invariant 
to the transformations of \g) by a time- dependent phase factor. This partial 
result, \g) {t) , 

\~g)(t)=e l Jo dt ' Wd ^\g) {t) (26) 
is a product between the periodic function \g)(t) and the phase factor 

exp(i f dt' (g\id t ,\g)) , 
Jo 

in general not periodic. Because the exact stationary solutions represent the 
limit case when \g) is time- independent, a periodic function \g)(t) is consid- 
ered as a good approximation for a constant one. This periodicity condition 
leads to a constraint on the phase [31]: 

j o dt'(g\td tl \g) = 2r m , T = — , (27) 

which gives for each n — 1,2, ... the amplitude u>o in eq. (22), undetermined 
up to now. 

Finally, the stationary wave function \Q) is obtained by projecting out 
the time independent Fourier component of \g)(t)' 

\ n ) = fL m ^ ■ (28) 

For applications to nuclei with a large number of particles, \g)(t) is chosen as 
a product of four factors: exp(— i4> p L v xp )\Z p )^ p ), exp(-i(f) n L v xn )\Z n ) ( CJn ), and 
exp(-i(p p L c xp )\a p ) (u)p) , exp(-i(f) n L xn )\a n ) {u)n) , corresponding to the valence 
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protons, valence neutrons, and core protons, core neutrons, respectively. Here 
\Z) are the HFB functions previously specified, and |a) = \ai, ct2, a Hc ) are 
the Slater determinants formed with n c single-particle states a^, A; = l,n c , 
occupied by the core nucleons. Both functions are periodic time-dependent 
through the oscillating frequencies oo p (t), w n (t). The expectation values re- 
quired in the quantization formula (27) are given by 

(Z T \Ll T \Z T )=2tr[iJ 2 (l+P)- 1 } , (29) 

(a T \L c XT \a T ) = 2tr[l x ] , (30) 

(Z T \d u \Z T ) = ±tr{[i,d u t\(l + t 2 )' 1 } , (31) 

(a T \d u \a T ) = , (32) 

where l x is the matrix of the single particle angular momentum operator l x . 

In the following only the A.V. state with n = 1 will be considered, ap- 
proximating the quantized phase of the function \g)(t) as Jq dt'(g\id t >\g) ~ fit. 
Thus, the state (28) becomes 

\n) = ±£dte™\g) {t) . (33) 

Its norm 

AT 2 = (m = ^ [ dt jT dt'e'^-^igitMt')) (34) 

will be calculated assuming small angle oscillation amplitudes a p , a n in eq. 
(20). In this case {g(t)\g(t')) has a simple expression 

(g(t)\g(t')} = Tl T=p , n w T (uj T (t) - u T (t'))W T (uj T (t) - u T (t')) , (35) 

where 

/ \ -lrj, , 7M ^M 1 +*( W )(^ ) t ] ,OC,\ 

w T (u = (Z T (u)Z T (0)) = ' ===== , (36) 

^det[l+t^]det[l + t&(t{01 

W T {u) = <a»|a T (0)> = [det(g-)t] 2 (37) 
i' = q° v %)q° v * , (38) 
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and q ov is the overlap matrix q°l = (* o (0)|* b (a;)) (A. 10). 

The state \ breaks the symmetries of the microscopic hamiltonian be- 
cause the projection of the particle and angular momentum quantum num- 
bers is approximate. In this form it can be used to evaluate only the total 
strength 



W\M(M1)\0)\- 



(39) 



of the magnetic transition from the + ground state |0) = \g(0)) to the A.V. 
mode uj r — 0, n — 1. The required matrix elements 



(n\MUMi)\o) = 

(n\Mis(Ml)\0) = 
for the isovector 



^<fdte- int (g(t)\M iv (Ml)\0) 
dte- iQt {g(t)\M is (Ml)\0) 



M iv (Ml) = 



g P - g n 




and the isoscalar 



MUM1) 



g P + g r , 




(L xp — L xn )/j,N 



->xp + L xn )nN 



(40) 
(41) 

(42) 
(43) 



components, respectively, of the transition operator Ai(Ml) = Ai is (Ml) + 
«Mj„(Ml) will be calculated by using the microscopic values g p — 1, g n — 0. 
In the same approximation of small amplitudes as above, we have: 



(g(t)\L xp \0) = w n (u n )W n (u n )[f p (u p ) + F p (u p )} , 
{g(t)\L xn \0) = w p (u p )W p (uj p )[f n (LJ n ) + F n (uj n )\ , 



(44) 
(45) 



where 



f T {u) = (Z T (u;)\Ll T \Z T (0)) = 2tr[%ii$(l+ii'J)- 1 } , 
F T (u) = (a T (u;)\L c XT \a T (0)) = 2W T (cj)tr[l x (q^)- 1 } 

(4)ob = (*oH|ix|*6(0)) • 

Because the matrices t are not known as analytical functions of a;, to avoid 
a numerical integration in (34), (40) and (41), the functions w, W, f, and 



(46) 
(47) 
(48) 



10 



F obtained numerically will be interpolated by polynomials in u, and in the 
end all the integrals will be calculated analytically. 

IV. Numerical Results 

a) The ground state band 

In order to obtain the energy levels of the ground state band in rare-earths 
nuclei it is necessary to solve the system of equations (9), (10) numerically. In 
practice, it turns out to be convenient to restrict the SO(2n) (HFB) treatment 
to a small number (10-12) of valence nucleons, while the other particles are 
considered in a HF state. Such a type of calculation was done for ^fGd, 
^fGd and 1 g|Er, considering 12, 12, 6 valence protons and 10,12,0 valence 
neutrons, distributed on 12 and 9 single-particle states, respectively. These 
sets of states were chosen to contain the Fermi level at uo = 0, and also to 
have no crossings between their energy levels and the levels from outside, 
when uj increases. The constants of the model u>i, u 2 , uj 3 (A.l) and G (3) 
can be expressed in terms of two parameters: the mass number A and the 
deformation 5. Using the standard formulas [35], Ui, i — 1,3 are separately 
given for protons and neutrons by the equations 

UJi = UJ 2 = U S {1 + -) , UJ 3 = U s {\ - ^5) , UJiUJ 2 ^3 = Ugr > ( 49 ) 

with u sp = 39.8A" 1 / 3 MeV and u sn = 44.8A- 1 / 3 MeV; 5 has the values 
determined by Lamm [35] 

5i56 Gd = 0.239 , 5i 5 s Gd = 0.250 , 5i68 Er = 0.273 , (50) 

to have an optimal fit of the single-particle energy levels, and G = 23 /A 
MeV. 

The calculated rotational energy levels (eq. (13), (14)) and the experi- 
mental ones are presented in figure 1. For 156 Gd, a set of slightly different 
parameters: 5 = 0.237, u sp = u sn = 41A -1 / 3 MeV, yield almost the same 
levels. However, a much different spectrum is obtained if the ground-state 
wave function is a Slater determinant of GO. single-particle states (the HF 
approximation). As it can be seen from figure 2, in a HFB state the in- 
dependence of the angular momentum for the valence nucleons of 156 Gd is 
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completely different than the one obtained by a HF calculation (the results 
for protons and neutrons in this case are practically the same). The HFB 
curves show a "diamagnetic" behavior [15] correlated with the variation of 
the order parameters A p , A n represented in figure 3. The existence of the 
superfluid layer also leads to differences between the kinematical and dynam- 
ical moments of inertia extracted by interpolation from the numerical results 
(with 5 = 0.239, u sp = 39.8A' 1 / 3 MeV, uj sn = 44.8A" 1 / 3 MeV) obtained for 
the expectation values L t {uj) = (g\L XT \g), r = p, n of the angular momenta: 

L p (uj) = 19.8cj + 2(k> 2 , L n (uj) = 33.7c; + 6.2a; 2 , (51) 

and of the energy 8 T (oo) = E T (uo) — E T (0), 

37 9 44 1 

Ep(u) = — - c; 2 + 0.31|c;| , £ n (u) = — -uo 2 + 0.54|oj| , (52) 

( [L T ] = h, [£ T ] = MeV). Such differences practically do not appear in a HF 
calculation, when: 

L p {uj) = 3AAu , L n (uj) = A3.2u , (53) 

^ / x 34.2 2 „ . , 43.2 2 . JN 

8^) = —^ , £ n (to) = —uj 2 . (54) 

It is interesting to remark that in all the cases investigated, the solution t 
of the system (9), (10) has the same symmetry as the matrix q, in the sense 
that q a b = implies t ab = 0. This feature strongly reduces the number of 
unknowns and also suggests that only the non-zero elements of t are the 
important degrees of freedom in a further time-dependent description of the 
pairing vibrations in rotating nuclei [23], [36] [37]. 

b) The angle vibrations of 156 Gd 

A direct numerical calculation shows that {Z T \dJ\Z^), r — p,n, are very 
small, and consequently the formula (27) becomes simpler: 

/ dt[(j) p L p (u p ) + (f) n L n (uj n )] = 2ir . (55) 

J 
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After the replacement of T and L T with their explicit formulas (21) (I p , I n 

are given by (52)) and (51), one obtains c; = ^2^/12.93, and the effective 
energy of the A.V. state 

E A . V . = = 1.57S1 . (56) 

The interpolation polynomials for the functions (36), (37), (46), (47) are 
explicitly 

w p (u) = 1 - 6.2c; 2 , W p (u) = 1 - 3.47c; 2 , (57) 

w n (u) = 1 - 8.7c; 2 + 27c/ , W n (u) = 1 - 3.6c; 2 , (58) 

f p (u) = w(16.97 - 9.64c;) , F p {u) = ujW p (uj)(7.QQ - 7.86c;) , (59) 

/ n (w) = 16.2w , F n (u) =u W n (u) (8.52 -9.5a;) . (60) 
The final results 

A/" 2 (c; ) = 2.7c; 2 - 17.56^ + 60.65^ - 112.22^ + 85.95c; 10 , (61) 

(n\M iv (Ml)\0) = (2.12w - 4.73^ + 4.0c^ - 1.34c^ + 0.13co 9 )fi N , (62) 

(tt\M is (Ml)\0) = (0.15c; - 0.26^ + 0.92c^ - 0.84^ + 0.13u%)fx N , (63) 

contain the constant u>o, undetermined only by the parameters 5 and A ap- 
pearing in the description of the ground state rotational band. In order to 
obtain its value, the fit of one characteristic quantity for the A.V. states with 
the experimental data becomes necessary. An interesting comment is that 
even before this fit, the B(Ml) iv and B(Ml) is strengths can be estimated in 
the limit of the small amplitude uo — > 0: 

B(Ml) iv \ U0 ^ = l.m f j? N , B(Ml) i8 \ U0 ^ = 0.009/4 • (64) 

Naturally the quantity which should be fitted is the elastic constant C, but 
as its value is not precisely known, the fit will be on the energy Ea.v., by 
using the experimental value of the dominant line from the spectrum, 3.1 
MeV. Then, all the other quantities, including C can be calculated: 

B(Ml) iv = 1.87/4 , B(Ml) is = 0.03/4 (65) 

a n = 7.33 deg , a p = 8.57 deg (66) 
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C = 79 MeV , u = 0.55 MeV 



(67) 



As it was expected [19] the value of C given above is about four times smaller 
than the value 311 MeV provided by the formulas of Palumbo [38]. The state 
\fl) proves to be excited from the ground state mainly by the isovector com- 
ponent of the M(M1) operator, with a strength close to the experimental 
value for the total strength of the 1 + fragments, B{Ml) ivfixp = 2.3 ± 0.5 n 2 N 
[38] [39]. The other results justify the approximation of small oscillation am- 
plitudes a p , a n . 

V. Summary and Conclusions 

The starting point in this work was the treatment of the pairing interac- 
tion in deformed rotating nuclei. Assuming that the pairing term is separable 
in the pp direction, and is determined only by a single-particle operator Q, 
independently on basis, the expression (1) was obtained. In fact, with Q = K 
this should also be found by rewriting the operator defined as usual in 
the spherical or Nilsson basis, if the sums could be made infinite. Arguments 
supporting this choice are: the good agreement with the experimental data 
obtained for the energy levels of the ground-state rotational band, and the 
ease of the calculation, as almost all quantities are known in closed analytical 
form. In applications the only quantity which must be obtained numerically, 
the matrix t, was proved to have most of its elements 0, as they are in the 
matrix q. This means in the case of a 12 x 12 matrix t that the number of 
unknowns (t ab , A) in the equation system (9), (10) is reduced from 79 to only 
18. Thus, by using the CO. approximation for the single-particle problem, 
by defining the pairing interaction using the time-reversal operator, and the 
SO(2n) coherent states approach to the many-body problem, we have an al- 
most analytical model for the treatment of the rotating superfluid systems. 
As a byproduct, the overlap coefficients q ov (A. 10) whose closed formulas 
were not found in the literature, have their own importance in connection 
to the representation theory of the Sp(4,R) group in the space of the 2- 
dimensional anisotropic harmonic oscillator states. Another result was the 
explicit form of the microscopic wave function for the A.V. states, including 
the pairing correlations. This was derived within the cranking model, in the 
adiabatic approximation for the motion of the proton and neutron deformed 
potentials. In the quantization procedure, beside the standard results on 
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the energy spectrum, it was also touched the more delicate problem of the 
stationary states. The prediction of the B(M1) strengths is in a good agree- 
ment with the experimental data, but of a greater importance appears to be 
the result on the characteristic of the model represented by the limit value of 
B(Ml) iv (64) which is in the range of the experimental values for rare-earths 
nuclei. This is an argument towards the physical relevance of the model in 
the A.V. description, making it interesting for further studies. 

Appendix: The cranked oscillator 

In the cranking model for the anisotropic harmonic oscillator potential, 
the single-particle part of the hamiltonian has the form 

3 , i 

h co = huk{a\a k + -) - ul x (A.l) 
k=i 1 

where 



4 = f£(*« - (A.2) 
V la muJk 

l x = ih[S(a,2as — a2a\) + D{a\a\ — 0203)] (A3) 

s = - w * + u * , D = ^-^ {AA) 

This can be diagonalized exactly by a canonical transformation to new boson 
operators A\, A k , 

A\ = a[ , A{= x ika) + yjkdj (A.5) 
3=2,3 

which satisfy, beside the canonical commutation relations, the equation 
[h co ,Al}=hQ k Al (A.6) 

W 2+ W 3 2 



Q 1 = u 1 , Q| 3 = i + uj 2 ± -yj(u% - u>l)* + ^{ul + ul) {A.l) 

It can be shown that in general, for the transformation (A.5) the scalar 
product (overlap) between the elements of the set of functions 

(A^) nk 

^ ni „ 2 „ 3 = njU ^=^000 , Afe^ooo = (A8) 
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with elements from the set 

$ nin2 „ 3 = n| =1 ^4 T $ 000 , a fe $ooo = (A9) 

both obtained by acting with the operators A\, and a\ on the vacuum states 
^ooo an d $ooo; respectively, is given by 

Qm 1 m,2rn3,n 1 n2n 3 = m 1 rn2in 3 \^ nm2n 3 ) (A. 10) 

= S mini y n 2 } .n 3 \m 2 \m 3 \ ^ £ X! 

k 2 =0 k 3 =0ji+j 2 +h +h=k 2 j[+j' 2 +l' 1 +l' 2 =k 3 


2k 2 n 2 p3fc 3 n 3 H rn2 + 'h+j' l -j2-j2,m 3 +h+V 1 -l2-l'2 



Xn h32hl2 n j' 1 j l 2l' 1 l l 2 



\j{m 2 + ji + J[ - n ~ Wijns + h + l[ ~k~ P 



The coefficient ?° 2n3 = ($on 2 n 3 |^ r ooo) has non-zero values only if n 2 , n 3 are 
both even or odd, and can be calculated easily in a harmonic oscillator rep- 
resentation for the operators a\, a k - If a k = + d/d£k)/yj2 then: 



$o„ 2 „ 3 = n fc=2 , 3 ^=g nfc fe) (a. ii) 



where 



*ooo = (^VV^-^-**. (A12) 

7T 

Mfc = 1/[1 + 2^(|y fcj | 2 - x,yyy] (A13) 
i=2 



(^23 - 2/23) (^32 - 2/32) - (x 22 - y 22 )(x 33 - y 33 ) 



and 



, n 2 \n 3 \^T 3 - (-l)^(Sf)^(T^) fc 



x t ^n-^~ m) iA15) 
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with 

„2 



T l-(-l) fc s 2 

h = o > ^ = — ; 



2 ' ' W 2 +l 

The factors R^ 2 i^ 2 are defined by the expansion 

n 

(4) B = »!E E ^^^(^^(aJ) 13 (^16) 
fc=0ii+i 2 +h+/2=fc 

and are, explicitly 

*& Ia = E (i-^)(i-^s)(-i)^ 

2(™- fc )/ 2 j 1 ! j2 !(f)!(|)!/ 1 !/ 2 ! 1 ' J 

The time- reversal operator K acts on $ nin2n3 through complex conjugation, 
and to obtain the matrix q, defined as qm 1 m 2 m :i ,n 1 n 2 n :i = {^m im2 m :i \K\^ nin2ns ) 
it is enough to use the expressions derived above for the case when (Aji)*, 
(Aj)* are written in terms of Ap Ay Because in this case 

(4)* = A\ , (4)* = £ a jk A) + p jk A, (A.18) 

3=2,3 

by choosing the phases such that x* k = (—l) j Xj k , y* k = (—l) j yj k , j = 2, 3, 
a jk = x 2 jX 2 k - y2jV2k + x 3k x 3j - y 3j y 3k {A. 19) 

Pjk = x 2j y* 2k - x* 2k y 2j + x 3j y* 3k - x* 3k y 3j (A20) 

are real numbers, we have s = 0, and ?° 2na are also real, non-zero only for 
n 2 , n 3 even, having the simple expression 



o 1 A-u 2 ] ^ ( l-u 3 ^ n^nsU/u^ 

<ln 2 n 3 (my^y} 1+U2 ) ^ 1 + U3 > ^2^+"3-2(l +U2 )(l +U3 ) • 1 • > 

It is interesting to remark that excepting for the term huji{a\ai + |) the 
CO. hamiltonian is a linear combination of the representation operators 
for the sp(4,R) Lie algebra {pjP^PjQk + qkPj,qjqk} j,k = 2,3 [40] [41] in 
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the Hilbert space Ti of the anisotropic harmonic oscillator states, and the 
transformation (A. 5) corresponds to a unitary representation operator U of 
the group Sp(4,R) in H: 

A\{u)=U{u)alU-\u) (A22) 

* ni « 2 » 3 (w) = U(uj)<S> nin2n3 (A23) 

whose matrix elements are the coefficients q ov . Moreover, the operators (A. 2) 
are connected to the boson operators for the isotropic harmonic oscillator 
b\ = ^JmuJ s /2h(xk — ipk/moOg) by the transformation a\ = V^^V^ 1 ^). 

Here V(6k) is an unitary operator generated by the element {{b\) 2 — b\)/2 of 
the representation for sp(2,R) Lie algebra in the Hilbert space of the harmonic 
oscillator states: 

V{B h ) = e e ^ 2 - b ^ 2 (A24) 

a\ = cosh 9 k b k — sinh# fe 5 fe (A25) 

cosh^ = i(/^+ & (A26) 

Denoting by 0o(£fc) the vacuum state for the b k bosons, and by <f) nk = 
KD~ 1/2 ( & !) nfe< M£fc), the "deformed" vacuum $000 becomes 

, ~ ,/(2n fc )!(tanh^) nfe 
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Figure Captions 

Fig. 1. The experimental and calculated energy levels of the ground state 
band for 156 Gd, 158 Gd, 168 Er. For 156 Gd it is shown also the HF result, 
neglecting the pairing interactions. 

Fig. 2. The angular momentum as a function of the rotation frequency for 
the considered valence nucleons of 156 Gd (5 = 0.237, cu sp = uj sn = 41A~ 1 / 3 
MeV). 

Fig. 3. The order parameters A p , A n for the valence nucleons of 156 Gd 
(5 = 0.237, cu sp = u sn = 41A~ 1/3 MeV). 
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